Robustness of non-abelian holonomic quantum gates against parametric noise 
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We present a numerical study of the robustness of a specific class of non-abelian holonomic 
quantum gates . We take into account the parametric noise due to stochastic fluctuations of the 
control fields which drive the time-dependent Hamiltonian along an adiabatic loop. The performance 
estimator used is the state fidelity between noiseless and noisy holonomic gates. We carry over our 
analysis with different correlation times and we find out that noisy holonomic gates seem to be close 
to the noiseless ones for 'short' and 'long' noise correlation times. This result can be interpreted 
as a consequence of to the geometric nature of the holonomic operator. Our simulation have been 
performed by using parameters relevant to the excitonic proposal for implementation of holonomic 
quantum computation [P. Solinas et al. Phys. Rev. B 67, 121307 (2003)] 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 

The use of uniquely quantum phenomena to process 
information has led to surprising results in quantum key 
distributions ID, information transfer protocols Q and 
computation [3j . From the point of view of the actual im- 
plementation of these theoretical protocols a main chal- 
lenge is posed by the fact that generically quantum states 
are very delicate objects quite difficult to control with the 
required accuracy. The interaction with the many degree 
of freedom of the environment causes the loss of informa- 
tion (decoherence) and moreover errors in processing the 
information lead to a wrong output state (control errors). 
The first problem has been extensively studied over the 
past few years and a few ways to overcome it have been 
proposed and experimentally realized. These strategics 
include error avoiding error correcting strategies Q 
and decoupling techniques . 

A new approach called topological quantum computa- 
tion has been argued to be able to effectively solve both of 
them and open new ways to inherently robust quantum 
computation 0. Information is encoded in topological 
degrees of freedom of a system which are not sensitive to 
the local environment-noise effects and then are robust 
against decoherence ||. Unfortunately to the date, no 
simple feasible physical system has been identified to this 
aim; in fact, the systems proposed are usually compli- 
cated many-particle ones living over a macroscopic non- 
trivial structure (e.g. torus or cylinder topology). On the 
other hand, we can develop topological information pro- 
cessing, where the operator used depends on topological 
controls and then are robust against a the unwanted fluc- 
tuations of the driving fields. In this case, a important in- 
termediate step is the geometrical quantum computation 
and particularly promising is the fully-geometrical ap- 
proach called Holonomic Quantum Computation (HQC) 

At variance with topological information processing, 



for geometrical QC several implementation proposals 
have been made; indeed the holonomic structure shows 
up in a variety of quantum systems, both in its Abelian 
(Berry) E3 and non- Abelian version El El El EI 
For this reason interest around HQC has recently growth 
leading also to proposals in which the adiabatic request 
can be relaxed and for implementing geometrical non- 
adiabatic quantum computation |15|. One of the appar- 
ent advantages of HQC is that the gating time for holo- 
nomic gates does not depends on the logical operator ap- 
plied but only on the adiabatic request; then HQC may 
lead to a new approach to implement complex operators 
difficult to construct with the standard dynamical gates 
(as discussed in E3)- ^ n view of its geometric nature 
i.e., dependence on areas spanned by loops in the control 
parameter manifold, HQC has been suggested to be ro- 
bust against some class of errors ^| E3 ■ Nevertheless 
thorough studies aimed to address this important issue 
are still relatively few and certainly not exhaustive [Isj . 

In this paper we will deal with the noise due to impre- 
cise control of the system parameters during the evolu- 
tion. This error source will be referred to as parametric 
noise. We will show that the fidelity of the holonomic 
gates displays three regimes, depending on the ratio be- 
tween the adiabatic time and the noise correlation time. 
This results can be understood in view of the geometri- 
cal dependence of the holonomic operator. We will study 
in detail the class of holonomic gate proposed in EH m 
which the physical system are semiconductor quantum 
dots, the logical qubits are excitonic quantum state con- 
trolled by ultrafast lasers. 

In section [n] after a brief review of the holonomic ap- 
proach, we describe the system used and the logical gate 
studied; moreover, it is discussed how we model the noise 
in the control parameters. In Sec. IIIII we give a descrip- 
tion of our simulations and show the results with different 
kind of noise processes for two single qubit gates and for a 
two qubit gate. A comparison with dynamical gates sub- 
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ject to the same noise is given too. Section llVl contains 
the conclusions. 



II. HOLONOMIC QUANTUM GATES WITH 
PARAMETRIC NOISE 

Let us consider a family T of isodegenerate Hamil- 
tonians H(X) depending on m dynamically controllable 
parameters A, in the HQC approach || one encodes 
the information in a n— fold degenerate eigenspace £ 
of an Hamiltonian H(Xq). Varying the A's and driv- 
ing adiabatically H(X) along a loop in the A manifold 
we produce a non-trivial transformation of the initial 
state \ipo) — ► U\i/jq). These transformations, known as 
holonomies, are the generalization of Berry's phase and 
can be c omp uted in terms of the Wilczek-Zee gauge con- 
nection [Tjj: U(C) = ~Pexp(§ c A) where C is the loop 
in the parameter space and A = A^dX^ is the 

u(n)— valued connection. If |Dj(A)) (i — l,...,n) are 
the instantaneous eigenstates of H(X), the connection 
is (AJ Q/3 = (D a \d/dW\D ) (a, (3 = l,...,n). The 
set of holonomies associate with a given connection is 
known to be a subgroup of the group of all possible n- 
dimensional unitary transformations; when the dimen- 
sion of this holonomy group and coincides with the di- 
mension of U (n) one is able to perform universal quan- 
tum computation with holonomies 0. 

For concreteness in this paper we will focus on the class 
of holonomic quantum gates analyzed in [l3l[T3 |. Logical 
qubits are given by polarized excitonic states controlled 
by femtosecond laser pulses. The parameters we have 
used in performing our simulations are those relevant to 
this specific kind of physical systems. 

First we concentrate to one-qubit gates. Despite this 
might look at first as major limitation, we observe that, 
as far as the holonomic structure is concerned, the two- 
qubit gates are very similar. So we expect that most of 
the results we are going to present here e.g., the existence 
of separate regimes, should, to a large extent, hold true 
for two-qubit gates too. 

The time-dependent interaction Hamiltonian in the in- 
teraction picture is 

H mt = -h(Cl + \E+) + Sl-\E£) + Slo\E%))(G\ 

+ h.c. (1) 

where \E\) (i = +, — , 0) are the polarized excitonic states 
(two logical and one ancilla) and \G) is the ground state 
(absence of exciton). This Hamiltonian family admits 
two dark states i.e., Hi nt (fl)\Di(il)) = 0, (i = 1,2). This 
two-fold degenerate manifold represent contains our en- 
coded logical qubit: \0) L := \E+), \1) L := \El). In Refs. 
[T2I IT^ ] it has been shown that the Wilczek-Zee connec- 
tion associated to the Hamiltonian family (JTJ allows to 
construct universal one-qubit gates. These are realized 
by giving an explicit prescription for driving the control 
parameter J7's along suitable adiabatic loops. 



We suppose now to add to the control Sl-field a 'small' 
noise which perturbs our trajectory on the control man- 
ifold. We test the robustness of the geometrical mix- 
ing single-qubit gate proposed in Ref. |T^ . To obtain 
this gate, we made the following loop in the parameter 
space : fi_(i) = £1 sinflcostp, = ^ sin 9 sin ip and 

Oo(t) = f2 cos 9 (where the il parameters is fixed and 
the 9(t) and tp(t) are time dependent) and the holonomic 
operator obtained at the end of the loop is U = e^" 7 " , 
(where ia y = \E^){Ej;\ - \El){El\). The geometrical 
parameter 4> — § sm 9d9dip is the solid angle spanned 
by the parameter vector f2 = (f2+, f2_, fio) on the pa- 
rameter manifold (sphere). Changing the relation be- 
tween 9 and tp we change the loop and then the value 
of <fi. We choose a loop in order to obtain <f> = it/ 2 and 
U - exp(ma y /2) - \E+){E£\ - \E-){E+\. 

The logical operator U depends only on geometrical 
parameters (i.e. solid angle swept on the parameters 
manifold), every perturbation that changes the trajec- 
tory in the control manifold changes the operator lead- 
ing to computational errors. This notwithstanding, the 
perturbations leaving (almost) unchanged this angle will 
not affect the holonomic operator. Then we expect that 
even strong fluctuations - provided their time scale is suf- 
ficiently fast - average out, leaving in this way the angle 
unchanged and thus not affecting the computation. 

The control parameters are the intensities and the 
phases of the lasers (since we suppose to be in resonant 
condition, we have fixed frequencies) perturbed by exter- 
nal noise (one for the phases and one for the intensities). 
To clarify which of these mainly affects the gate opera- 
tion, first we separate the two types of errors (intensity 
and phase fluctuations) and then we apply both of them. 

We use a straightforward model for the noise : we ex- 
tract a random number from a probability distribution, 
we add a constant noisy field for the time T n to the evolv- 
ing control field, then we extract another random number 
and so on. To simplify the simulations we choose T n in 
such a way that T a d is a multiple of T n (T a d = n r T n ). 
The fundamental parameter is the noise time T n that 
is the lapse of time of each random extraction; i.e. it 
represents the time scale of each random fluctuation. 

For the 'intensity' noise we modify only the value of 
the Rabi frequencies fi. We have three lasers turned on 
and we suppose they have independent fluctuations 5tti 
(i.e. every T„ we extract three random numbers). The 
evolution on the control manifold is described by: 

f2_(t) = Q sin cos ip + SO,- (t) 
fi + (t) = ^ sin 9 sin ip + 5fl + (t) 
Sl (t) = fl cos9 + Sfl (t) 

In this case, the Rabi frequencies remain real parame- 
ters, but, if we introduce a 'phase' noise (fl, — > e^flj), 
they acquire an imaginary part. The random numbers 
are taken as before and the evolution in the control man- 
ifold is: 
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FIG. 1: Fidelity for holonomic gate with — 50 fs, T ac i — 
7.5 ps (fi T ad = 150) and SO, = 0.1JX (A) Slowly varying 
fluctuations and (B) fast varying fluctuations. 
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FIG. 2: Loop in the parameter space for the holonomic 
gate with (solid line) and without noise (dashed line). With 
extraction of 2 (T n = T a d/2) (A) and 70 random numbers 
(Tn = Tad/70) (B). The sphere radius is normalized to Q = 1. 



fi_(t) = e 1 ^- £1 smOco&ip 
Q+(t) = £1 sin sin 
ft (t) = e iio ft cos 9 

The most general and complicated situation is when 
both 'intensity' and 'phase' noise are present. 

III. SIMULATIONS 

For all the simulation we choose the parameters used 
in Refs. 0, Q] which satisfy the adiabatic condition 
: ft = 0.02 fs" 1 , T ad = 7.5 ps and ft T ad = 150. The 
probability distribution for the noise is a Gaussian with 
zero mean and < a >=< ^ >= 0.1 where 5fl is the 
fluctuation of the Rabi frequency. Though this value is 
far below the experimental control, we use it to under- 
stand the robustness of holonomic gates against strong 
perturbations. 

Once the evolution of the state with the noise has 
been computed, we must compare it with the ideal one 
in which the noise is absent. In order to make such a 
comparison quantitative we introduce the fidelity 

r = V '(r tt\pout\r d ut ) (2) 



where iV'ott) ^ s * ne nna ^ state without noise and p out 
is the density matrix of associated to the noisy final 
state. In our case the evolution is unitary, then p out = 
\ipout se ) (' l Pout se \ an( i t ne fidelity reduce to be a scalar 
product between the noisy and the ideal state. To elimi- 
nate the dependence of on the initial states we make 
a sampling of the initial state space and then average 
the results. Even if we have a four dimensional working 
space, the initial state space has dimension two; in fact, 
in the ideal gates (once satisfied the adiabatic condition) 
we always start and end in a superposition of the logical 
states \Et) — \E^). This simplifies the sampling proce- 
dure because we can take the Bloch sphere as samp ling 
space. We sample the Bloch sphere with 18 states [2(j 
and for each of them calculate the fidelity. Fixed the ran- 
dom numbers extracted during the evolution (i.e. fixed 
T„), we have also to take into account the dependence 
from the series of the random number extracted. For ev- 
ery sampled state we make five different realization with 
different extracted random numbers and calculate the av- 
erage fidelity for the particular initial state. Finally, we 
make the average fidelity for all the sampled states to 
obtain the fidelity of our gate for a fixed T n . 

In the Figures presented here on the y axis is plot- 
ted the fidelity and on the x axis is plotted the number 
of extractions (that is the the ratio T a d/T n between the 
adiabatic and noise time). 
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FIG. 3: Loop in the parameter space for the holonomic gate 
with noise and the extraction of 5000 random numbers (T n = 
Tad/5000). The sphere radius is normalized to fi = 1. 



In Figures ^ w e report the simulated evolution for two 
regimes of noise. In fig. Q (A) we plot the fidelity when 
we extract up to 100 random numbers during the evolu- 
tion (T a d/ 100 < T n < T a d)- Up to 30 random numbers 
the average fidelity is 0.875, while it decreases up to a 
minimum of about 0.5 (with total average of 0.632) if we 
extract more random numbers. 

A possible interpretation of this effect can be given 
looking at Figure [21 where we show the evolution on the 
parameter sphere. In figure |2| (A) n r — 2 (T n = T a d/2) 
and we change the noise field twice during the evolution. 
Despite the intense noise, the shape of the loop is still 
clearly visible, it is simply shifted with respect to the 
ideal one. The value of the resulting solid angle swept is 
near to the ideal one. In figure |2 (B) we extract n r — 
70 random extractions during the adiabatic evolution 
(T n = T a d/70). The fluctuations are too intense and too 
few to cancel out then, as we expected, the solid angles, 
swept, respectively during the ideal and the noisy loops, 
are different. This can explain the result in figure ^ (A). 

As stated before, even if we have strong fluctuations, 
we expect that if we extract many random numbers the 
noise in average does not affect the solid angle and leave 
the holonomic operator unaffected. This is what seems 
to be confirmed by the simulations illustrated in Figure^ 
(B) where we extract from 50 to 5000 random numbers. 
The fidelity increases and it is even better of those in fig- 
ure A); the average fidelity is 0.918 but it increases to 
0.956 if take into account the values from 1000 to 5000 
random extractions. The relative loop in the parameter 
space is shown in figure where we extract 5000 ran- 
dom numbers during the evolution: the fluctuations have 
<50 = 0.1 f2 but are so quick that cancel out. 

From these simulations it is evident that we have three 
regimes in our model which can be explained as the sign 
of the geometrical dependence of the holonomic operator: 

• Slowly varying random fluctuations (T a d/T n * 1): 
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FIG. 4: Fidelity for small perturbations 8Q = 0.01S1. Param- 
eters as in figure Q 

the loop basically maintains its shape and it is sim- 
ply shifted. This situation does not affect the gate 
too much. 

• Intermediate regime (50 < T a d/T n < 100): the 
intense fluctuations badly modify modify the loop 
shape and alter the gate operator. 

• Fast varying random fluctuations (T a d/T n 3> 1): 
the fluctuations effectively cancel out and do not 
change the operator. 

These geometrical features are independent from the 
ratio <5f2/f2 and persist even for small values of the fluc- 
tuation Sft. If we decrease the intensity of the noise our 
adiabatic gate improves as shown in Figures Q where, 
with the same parameters (i.e. Rabi frequency and adi- 
abatic time), we choose S£l = 0.01S1. We note that the 
fidelity increases but the features are the same of the pre- 
vious simulations (i.e. three regimes). 

In Figure (J5J the noise is applied only to the phase of 
the control field with the same adiabatic parameters of 
the previous simulations. Since the way of producing the 
operator is always the same (i.e. loop in the parameter 
space with a noise), we expect that the effects of noise 
are the same of those for 'intensity' noise. This can be 
clearly seen in figures JSJ where we find the same features 
of the previous plots. The fidelity is much better respect 
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to 'intensity' noise and we can say the the 'phase' noise 
does not affect our gate. 

We finally discuss the case of both intensity and phase 
noise (Figure EJ). As we expect, the main part of the 
error is given by the intensity noise . The three regimes, 
discussed previously, are evident also in this case. 

In Figures Q we show the population of the non- 
logical states (\G) and \E^)) at the end of the gate appli- 
cation as function of T ac t/T n . In an ideal adiabatic gate 
(for T a< i — ► oo ) these states are not populated. In Figure 
Q (A) these populations increase with the T a d/T n ratio 
due to the strong and fast fluctuations of the Rabi fre- 
quencies which perturb the Hamiltonian. For very fast 
varying fluctuations (Figure (B)) the undesired pop- 
ulations decrease because of cancellation effects. From 
this analysis we can say that for slowly and fast varying 
random fields (i.e. the two extreme regimes) we remain 
in the logical computational space. 

We wish to compare this holonomic gate with a stan- 
dard dynamical gate, the latter one being characterized 
by the same unitary operator of the holonomic one and 
having different (typically shorter) gating time (the time 
need for the application of the gate). In our system 
this means to consider as logical states \G) = |0) and 
\E l L ) = |1) and apply a 7r— pulse laser sequence to pro- 
duce a transition \G) <-» \E\). 

We made similar simulations for dynamical gate with 
the same parameter (fl = 0.02 fs _1 ) and with the same 



noise. 

In Figures JSJl we make a comparison between holo- 
nomic (solid line) and dynamical gates (dashed line) sub- 
ject to the same intensity and phase noise. The compari- 
son is not direct since the gating times are different. This 
means that the ratios T a d/T n and Td yn /T n are different 
for adiabatic and dynamical gates. To compare the ef- 
fect of the gate subject to the same noisy field (i.e. with 
the same noise time T n ), we have to take into account 
that the dynamical gates, in our model, are about 100 
times faster than the adiabatic ones (seeRefs. [H Q). 
This means that if during the adiabatic evolution the 
noise changes n^ d = T a d/T n times, for the dynamical 
gates it changes only nf yn — T a£ j/(100 T n ) times (if for 
the dynamical noise nf yn < 1, it changes only once). In 
other words, if for dynamical gate we extract n dyn ran- 
dom numbers during the evolution, for the holonomic 
gate we extract n% d = 100 n dyn random numbers. This 
means that to compare the fidelities we have to look in 
the fast varying fluctuating noise (T a d/T n 1) region of 
the previous figures. 

For this reason in Figures (A) and (B) we put two 
different scales for the adiabatic (T a d/T n ) and dynamical 
{Tdyn/Tn) gates and the plots have been rescaled taking 
into account for the different gating time. 

Numerical simulations in Fig. [SJshow that the perfor- 
mance of holonomic gates and dynamical gates are com- 
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FIG. 8: Comparison between holonomic (solid line) and dy- 
namical (dashed line) gates with (A) S£l = O.lfi and (B) with 
SQ — 0.01f2 . On the top T a d/T n for the holonomic gate and 
on the bottom Td yrL /T n for dynamical gates are reported. 



parable in the region where the first ones are reliable; 
that when T ad /T n ^> 1 (see discussion above). 

Both dynamical and holonomic gates can be further 
improved. Since the dynamical gates are not subject to 
adiabatic constraints, we can choose different parameters 
in order to minimize the effect of the noise but this can 
affect the gating time. 

For holonomic gates, given a noise with fixed correla- 
tion noise time T n , we can try to change adiabatic time 
in order to modify the T ac i/T n ratio. This should allow us 
to fall in a 'good' regime i.e. fast or slowly varying fluc- 
tuations. Decreasing adiabatic time to get in the small 
T a d/T n region can produce new errors due to the lack 
of adiabaticity during the evolution and then it must be 
treat carefully. Increasing the adiabatic time to enter in 
the region of great T a d/T n leads both to a better pre- 
cision in the adiabatic gates and to the cancellation of 
noisy fluctuations but gives longer gating times. 

To complete the set of universal quantum gates, we 
need another single qubit gate and a two qubit gate. 

In the first case, we apply an intensity noise to the 
phase gate presented in Ref. 0] with = 0, £1+ = 
-fisin(0/2) e lv and n = Qcos(0/2). The holonomic 
operator is U = e"^ B ^ I where <j> = \ § smOdOdip. 
The results are shown in Figures © and the structures 
discussed above are evident. 

Finally, in Figures l|l(J|l we present the simulations of 



the phase shift two qubit gate in Ref. |13l | with an in- 
tensity noise. In this case we use two exciton states 
{\E l L E J L ) and \GG)) and a two photon interaction Hamil- 
tonin similar in structure to Q but with Rabi frequency 
ft e ff = 2hil 2 /5 (where 5 is the laser detuning we need 
to avoid single photon transition and create two exciton 
states). 

In Figure l|10|) (A) where the second regime (where the 
fidelity decreases) is present but less evident and seems 
to be compressed in the slowly varying random fluctua- 
tions zone (10 < T a d/T n < 30); for great T a d/T n ratio 
the fidelity decisively increases as in the previous figures 
(Figure (|10|) (B)). Moreover, we note that the fidelity be- 
tween the adiabatic final states with and without noise is 
high and close to 1 even if we have chosen the adiabatic 
parameter (O e /r T) smaller than the one used for sin- 
gle qubit gates. These can be explained as consequence 
of the fact that in the imperfect adiabatic evolution un- 
wanted states gets populated. On the other hand the 
effect of the fluctuating noise is to induce undesired tran- 
sition as well. These two effects are superposed and with 
smaller adiabatic parameter the effect of the noise seems 
to be less important. 
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FIG. 9: Fidelity for the single qubit phase shift gate with to 
intensity noise. The parameter are the same of figure Q 

IV. CONCLUSIONS 

We numerically studied the robustness of a non-abelian 
holonomic quantum gate against stochastic errors in con- 
trol parameters. The robustness of logical gates show 
three regimes upon the variations of the noise correla- 
tion time T n . A possible interpretation of these regimes 
can be given on the basis of the the geometric (i.e. solid 
angle swept in the parameter space) dependence of the 
holonomic operator. For fast random varying fluctua- 
tions we have a good robustness of the holonomic gate 
since, as argued in other papers 0, 0, the fluctua- 
tions during the loop tend to cancel out. For random 
varying fluctuations in the intermediate regime the holo- 
nomic gates are significantly corrupted because the fluc- 
tuations strongly deform the parameter loop. For slowly 
random varying fluctuations the performance improves 
again. This fact is not surprising as it may seem, indeed 
the loop in the parameter space turns out in this case 
to be simply is shifted rather than deformed; then simi- 
lar solid angles are swept. Our analysis suggests that the 
main noise source is given by fluctuations in the intensity 
of the control parameters whereas the phase fluctuations 
do not seem to sizeably affect the gate studied. The effect 
of the noise decreases with the variance of the intensity of 



the fluctuations and for SQ/fl — 0.01 we have average fi- 
delity close to 1. A first comparison shows that holonomic 
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FIG. 10: Fidelity for the two qubit phase shift gate with to 
intensity noise. Parameters are : 8 = 5 meV, = 5/15 and 
T ad = 0.8 ns. 



and dynamical gates have comparable performance in the 
Tad/T n > 1 region. 

We performed similar simulations for two single qubit 
gates and for a two qubit gate in order to complete the set 
of universal quantum gates. For the single qubit gates we 
obtain similar results. For the two qubit gate the three 
regimes descibed above are present but less evident. 

We believe that our analysis and conclusions should 
be extended rather easily to different sort of system pro- 
posed for HQC: for example, it definitely extends to the 
model proposed in ref. [12 ] since the involved holonomic 
structure is isomorphic to the one here analyzed. The 
general features of our results should hold also more gen- 
eral situations, since they apparently do not rely on the 
detailed features of Hamiltonian , rather on the gen- 
eral structure of holonomic computations. A related, 
though logically distinct, issue is the robustness of HQC 
against environmental decoherence [2j]]. This is clearly 
an important subject to be addressed in future investi- 
gations. 

P.Z. gratefully acknowledges financial support by 
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Union project TOPQIP (Contract IST-2001-39215) 
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